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Abstract 

The Kramers problem in the energy-diffusion hmited regime of very low friction is 
difficult to deal with analytically because of the repeated recrossings of the barrier that 
typically occur before an asymptotic rate constant is achieved. Thus, the transmission 
coefficient of particles over the potential barrier undergoes oscillatory behavior in time 
before settling into a steady state. Recently, Kohen and Tannor [D. Kohen and D. J. 



Tannor, J. Chem. Phys. 103, 6013 (1995)] developed a method based on the phase 
space distribution function to calculate the transmission coefficient as a function of 
time in the high-friction regime. Here we formulate a parallel method for the low- 
friction regime. We find analytic results for the full time and temperature dependence 
of the rate coefficient in this regime. Our low-friction result at long times reproduces 
the equilibrium result of Kramers at very low friction and extends it to higher friction 
and lower temperatures below the turn-over region. Our results indicate that the single 
most important quantity in determining the entire time evolution of the transmission 
coefficient is the rate of energy loss of a particle that starts above the barrier. We test 
our results, as well as those of Kohen and Tannor for the Kramers problem, against 
detailed numerical simulations. 

1 Introduction 

The classic work of Kramers in 1940 |l| on reaction rates in which the effect of the solvent 
was taken into account in the form of a Markovian dissipation and Gaussian delta-correlated 
fluctuations has spawned an enormous and important literature that continues to flourish 
10, ^. Recent advances in experimental methods that make it possible to monitor the 
progress of chemical reactions on microscopic (even single-molecule single-event) spatial and 
temporal scales |^ have led to a revival of interest in theoretical approaches that examine 
the problem on these detailed scales. 

The enormous literature on the theoretical front since the appearance of Kramers' seminal 
paper has evolved in many directions that include more formal derivations of Kramers' own 
results, extensions to larger parameter regimes and to non-Markovian dissipation models, 
generalizations to more complex potentials and to many degrees of freedom, analysis of 
quantum effects, and apphcation to speciflc experimental systems. 

One such recent direction, developed by Kohen and Tannor (KT) [^, deals with the 
derivation of the rate coefficient in the Kramers problem and in the more general Grote- 
Hynes problem (that is, the Kramers problem extended to a non-Markovian dissipative 
memory kernel) so as to obtain not only the asymptotic rate constant but the behavior of 
the rate coefficient at all times. This derivation is based on the reactive flux method [@, |^, H 
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and allows the interpretation of the salient features of the time dependent rate coefficient 
k{t) in terms of the time evolution of the representative distribution functions that originate 
at the top of the barrier, thus paralleling closely and usefully the methods used in numerical 
simulations of the problem. KT analyze in detail the dependence on the dissipative memory 
kernel of the time it takes the rate coefficient to reach its stationary (equilibrium) value. 
Their extensive analysis, however, does not cover a number of parameter regimes, nor do 
they check their time-dependent results against numerical simulations. We undertake such 
studies to complement their work. Thus, we extend the parameter regime of analysis and 
check some of their results as well as our new ones against numerical simulations. 

In this paper we restrict ourselves to the Markovian regime, that is, to the original 
Kramers problem, deferring the analysis of the non-Markovian case to a subsequent pub- 
lication. The dynamical equation of interest to us is thus the generic reaction coordinate 
problem 

P 

<1 = — 

m 

dV{q) , . . . 

P=-1P 

where q{t) is the time-dependent reaction coordinate, a dot denotes a time derivative, 7 
is the dissipation parameter, V{q) is the potential energy, and fit) represents Gaussian 
delta-correlated fluctuations that satisfy the fluctuation-dissipation relation 

{f{t)f{r)) = 2^kBT5{t-T). (2) 

A;^ is Boltzmann's constant and T is the temperature. The potential V{q) is a double- well 
potential that is often (and here as well) taken to be of the form 

V(,)^V'o(^4 + i)^^tf-lf (3) 

(see Fig. |1|). The minima of the potential occur at g = ±1, the potential barrier has a 
maximum at g = 0, and the barrier height there is Vo/4. The parameter Vq can be used as 
the unit of energy, and henceforth we set it equal to unity. The barrier height is assumed to 
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be large compared to the temperature (i.e., fc^T ^ 1/4); otherwise the notion of a barrier 
crossing process loses its meaning. 



V(q) 




Figure 1: Potential energy for the Kramers problem as given in Eq. (^) with Vq = 1. The 
minima of the potential occur at g = ±1, the maximum at g = 0, and the barrier height is 
1/4. e measures the energy above the barrier, is the time it takes a particle of energy e 
above the barrier starting at g = to return to g = 0, and is the net energy loss in this 
half orbit. The quantities and n are discussed in Section |^. 

The problem of interest is the rate coefficient k{t) for an ensemble of particles whose 
positions evolve as realizations of g(t). The coefficient k{t) measures the mean rate of 
passage of the ensemble across g = from one well to the other. This crossing rate and, in 
particular, its asymptotic value k = k{oo), is associated with the rate constant of the process 
represented by the reaction coordinate. One usually focuses on the corrections to the rate 
obtained from transition state theory (TST) and therefore writes 

k{t) = (4) 

where k'^^^ is the rate obtained from transition state theory for activated crossing, which 
for our problem and in our units is 

TT 
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The deviations from this rate constant are then contained in the transmission coefficient 
K,{t). The construction of K,{t) is discussed in subsequent sections. In the remainder of this 
work we designate the stationary or asymptotic transmission coefficient K,{t —>■ oo) as Kgt- 




Figure 2: Transmission coefficient Kgt vs dissipation parameter 7 for two temperatures ob- 
tained from direct simulation of Eq. (0) as described in Section |. Sohd circles: ksT = 0.025; 
triangles: fc^T = 0.05. As already predicted by Kramers, the high friction result (short- 
dashed curve) is independent of temperature and equal to —(7/2) + ((7^/4) + 1)^^^, which 
approaches 7"^ as 7 ^ 00. At very low friction Kramers predicted that is proportional 
to 'j/ksT. The dashed and solid curves, based on our predictions for the behavior of Kst 
beyond the very earliest linear regime, are discussed in Section |^. 

The input parameters of the problem are the dissipation (7) and the temperature [ksT). 
The behavior of the stationary-state transmission coefficient function of these pa- 

rameters, particularly as a function of 7, was the subject of a great deal of work in the 80's. 
It was already realized by Kramers that the behavior of Ust vs 7 would exhibit a maximum 
(cf . Fig. 1^) since at intermediate and high dissipation ( "diffusion limited regime" ) the rate 
constant decreases with increasing dissipation while at low dissipation ("energy-diffusion 
limited regime") it increases with increasing dissipation. Indeed, Kramers found that at as 
7 — > CX3, ^ 1/7 while as 7 — 0, ~ j/ksT. These limits were obtained from separate 
arguments, however, which therefore gave no indication of where precisely the "turnover" 
(maximum) might occur. Subsequent work attempted to address this problem [^. The other 
parameter dependence, namely, on temperature, has also been the object of study. For ex- 
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ample, it is understood that the high-dissipation Kgt is independent of temperature whereas 
Kst does depend on temperature in the low-dissipation regime. This is clearly seen in the 
simulation results in Fig. 0. Discussion of Fig. |^ will be resumed later. 

The analysis of KT has carried the problem further in that they actually dealt with the 
full time-dependent evolution of the transmission coefficient. They derived an expression for 
K(t) in terms of the time-dependent phase space density, assumed a particular form for this 
density, and studied how goes to its steady state value as the phase space density relaxes 
to equilibrium. In carrying out this program, however, they relied on approximations that 
are appropriate for the diffusion limited regime and therefore were able to calculate K{t) only 
in this regime. They did not check their time- dependent results against simulations, although 
they did confirm numerically that the (temperature-independent) steady state values for the 
transmission coefficient that they predict in this regime are satisfactory. They also did not 
deal with the energy-diffusion limited regime. 

Thus, to complement the work of KT in the Markovian regime we accomplish three goals 
in this paper: 

1. We carry out simulations of the time-dependent transmission coefficient to assess and 
confirm the validity of the KT formulation in the diffusion-limited regime. Their results 
agree very well with simulations in this regime. 

2. We formulate the complementary theory for the time-dependent transmission coeffi- 
cient in the energy-diffusion limited regime and compare our results with numerical 
simulations in this regime. Our theory captures the complex oscillatory and stationary 
state behavior of K{t) in this regime very well. 

3. We calculate and check via our simulation results the temperature dependence of the 
asymptotic transmission coefficient in the energy-limited diffusion regime beyond the 
Kramers result. 

In Section ^ we briefly review the reactive flux formalism. Section ^ contains a short 
description of the numerical methods used in our simulations. The results of KT in the 
diffusion-limited regime are recalled and compared to our simulation results in Section ^ 
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The departure of the KT results from the simulations at low dissipation are illustrated in 
this section. In Sections ^ and |^ we present our theory for the time-dependent transmis- 
sion coefficient K{t) and its asymptotic value in the energy-diffusion limited regime. The 
comparison of our results with those of numerical simulations is presented in Section |^. We 
conclude in Section E. 



2 Reactive Flux Formalism 

Two decades ago saw the development of the reactive flux formalism for the rate constant 
over a barrier |^. This formalism was important because it made possible efficient 

numerical simulations of the rate constant without having to wait the inordinately long time 
that it takes a particle at the bottom of one well to climb up to the top of the barrier. In 
the reactive flux method the problem is formulated in terms of the particles of the thermal 
distribution that are sufficiently energetic to be above the barrier at the outset. In this way, 
it is not necessary to wait for low-energy particles (the vast majority of all the particles) to 
first acquire sufficiently high energies via thermal fluctuations. 

We follow the notation of KT. The rate constant k{t) in the reactive flux formalism is 



6p{qo)ep[q{qo,Vo,t] 

_ {vo5{qo)dp[q{qo,Vo,t)]) 
iORiqo)) 



(6) 



where the top of the barrier is at position go = 0, Op{x) = 1 if x < and otherwise, and 
Op = 1 — 6 p. The brackets ( ) represents an average over initial equilibrium conditions. With 
the particular choices made in Eq. one is calculating the transition rate from the left well 
to the right. KT proceed through a series of steps that flnally yield for the transmission 
coefficient introduced in Eq. (11) the relation 
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where 



dq dvW{q,v,t] qo = 0,Vo) 



(8) 



and W{q, v, t; go = 0, Vq) is the conditional phase space distribution function that corresponds 
to an ensemble of particles starting at (go = 0, Vq) at t = 0. 

The (nonequilibrium) conditional probability distribution W clearly lies at the crux of 
the calculation: it is the distribution associated with a Langevin equation or a generalized 
Langevin equation that describes the evolution of the position q{t) and momentum p{t) = 
mv{t) of the ensemble of particles. In the problem of interest here, W is the distribution 
associated with Eq. (|l]), that is, it is the solution of the Fokker-Planck equation (we set the 
mass equal to unity) 



dW _ _d_ ^x^^ 
dt dq dp 



+ 'ykBT—-^. 



(9) 



with appropriate initial and boundary conditions. 



3 Simulations 

The numerical solution of Eqs. (P is performed according to the following main steps. The 
integration is carried out using the second order Heun's algorithm 0, which has been tested 
in different stochastic problems with very reliable results |TD|. A very small time step is used. 



ranging from 0.001 to 0.0001, as in Ref. 0. The numerical evaluation of the transmission 
coefficient nit) follows the description of Refs. 0, |^. We start the simulation with particles 
(either 4000 or 1000 depending of the circumstances), all of them above the barrier at g = 0, 
half with a positive velocity distributed according to the Boltzmann distribution in energy, 
which translates to the velocity distribution f exp(— f ^/2A;bT), and the other half with the 
same distribution but with negative velocities. 

The transmission coefficient is extracted from these sets of simulated data by calculating 



8 



N+{0) A^_(0)' 



(10) 



where A^+(t) and N^{t) are the particles that started with positive velocities and negative 
velocities respectively and at time t are in or over the right hand well (i.e. the particles for 
which q{t) > 0). 




Figure 3: Transmission coefficient K,{t) vs time for two values of the dissipation parameter 
7 and for temperature fc^T = 0.025. The solid curves are the results of our numerical 
simulations. For high dissipation (7 = 1.0, on the right of the maximum in Fig. 0), the 
transmission coefficient decreases monotonically. For low dissipation (7 = 0.05, on the left 
of the maximum in Fig. |^) the decay is oscillatory. The dashed curves are those calculated 
from the theory of KT (cf. Eq. ([I^), which is seen to predict the numerical results very well 
for high dissipation, the regime for which the theory is designed. The agreement seen in the 
figure is typical for large dissipation parameters. Not surprisingly, the theory fails at low 
dissipation. 

Figure |^ exhibits examples of the temporal behavior of the transmission coefficient, 
Eq. ([I0|), for two typical different parameter regimes as described in the figure caption. 
Of interest in this paper is the low-dissipation oscillatory behavior, which is illustrated and 
discussed further in subsequent sections. 

As discussed by Straub et al. whereas the exact transmission coefficient reaches a 
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constant non-zero value that corresponds to the equihbrium transition rate for the problem, 
the transmission coefficient calculated using the reactive flux method flattens out but con- 
tinues to decrease with time. If the temperature is low, this decrease is slow and the results 
appear flat on the sort of time scale shown in Fig. 0. In this case the value of K,{t) in this 
flat region is identified with Kgt- If the temperature is not so low, then the decaying tail is 
extrapolated back to its intersection with the vertical axis according to the relation 

n{t) ~ Kste-""' (11) 



where K is the decay constant of the tail. Thus, the value of the intersection is identified as 

The numerical values of Hgt presented in Fig. ^ are obtained according to these schemes. 
Numerical uncertainties are fairly large and difficult to avoid for small values of Kst- 



4 Diffusion Limited Regime 

Equation (|^) can not be solved exactly for the bistable potential, and therefore the dis- 
tribution W is the quantity in the calculation of the transmission coefficient that requires 
approximation. The quality of the result of course depends on the quality of the approxima- 
tions. KT focus on the regime of moderate to high dissipation, that is, on the Grote-Hynes 
model. In this regime it is appropriate to adapt the theory of Adelman |ll| to the assumption 



that the barrier is parabolic. With this approximation for the conditional probability, KT 
derive expressions for n{t) from which the Grote-Hynes results are recovered in the long-time 
limit. Their result is 

<t) = ; ^^2) 



where 



I ^2 ry 
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and uj is the "frequency" of the top of the barrier, which for the potential is unity. 

In order to assess the KT results and their regime of validity, we point again to Fig. ^ 
where we show the asymptotic transmission coefficient as a function of the dissipation param- 
eter obtained from simulations based on Eq. (|^. In this figure we have exhibited the results 
for two temperatures. The high dissipation side of the curve is independent of temperature, 
but the low dissipation transmission coefficient is inversely proportional to temperature. The 
short-dashed curve is the Kramers/Grote-Hynes result (and therefore also the result recov- 
ered by KT). Clearly, it is excellent at high dissipation (a well-known result), but does not 
capture the turnover and subsequent decrease in the transmission coefficient as one enters 
the energy-diffusion limited regime. 

Figure |^ exhibits a typical simulation result in the high dissipation regime (7 = 1.0) 
(monotonic somewhat jagged curve) and also the associated KT prediction Eq. ([T^) (dashed 
curve). The agreement is seen to be very good over the full time-dependence of the trans- 
mission coefficient. The very small discrepancy in the asymptotic regime is also typical, with 
the theoretical results slightly above those of the simulation. This may in part be due to the 
issues surrounding Eq. (0). In any case, this degree of agreement between the KT theory 
and the simulations is typical in the regime well beyond the turnover in Fig. |^. The theory 
fails, however, at low dissipation parameters. In Fig. |] we show typical results in this regime 
- these are for 7 = 0.05 and ksT = 0.025. Although the two values of 7 that we have chosen 
for this figure lead to similar values of the stationary transmission coefficient, the decay 
towards these values is clearly very different in the two cases. The simulation results (con- 
tinuous curve) and the KT result (dotted curve) agree at early times, until the simulation 
results rather suddenly drop and settle in an oscillatory way to a completely different lower 
asymptotic value than the theoretical result [0. The KT result in fact asymptotes to the 
Grote-Hynes rate, which in this regime is not correct. In the next sections we obtain results 
that predict the correct oscillatory and asymptotic behavior of the transmission coefficient 
in the very low dissipation regime. 
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5 Energy-Diffusion Limited Regime - Orbits 



In this section and the next, we detail the approach complementary to that of KT to obtain 
analytic results for the time-dependent transmission coefficient in the energy- diffusion limited 
regime. Our results are compared to numerical simulations in Section |^. 

For our calculation it is convenient to write the transmission coefficient as (cf. Eq. (|10D) 

«:(t) = «:+(t)-«:_(t), (14) 

where K,+ {t) is the fraction of particles that started above the barrier with positive velocity 
f > whose coordinate q{t) > 0, i.e., those that are over or in the right hand well at time 
t given that they were over the right hand well initially, and k_ {t) the fraction that started 
above the barrier with negative velocity f o < and are over or in the right hand well at time 
t. 

When dissipation is very slow, the particles that start above the barrier lose their energy 
very slowly as they orbit around at an almost constant energy during each orbit (recrossing 
the barrier many times if they start with sufficiently high energy). One can calculate the 
approximate energy loss per orbit and thus keep track of how long it takes a particle of a 
given initial energy above the barrier to lose enough energy to be caught in one well or the 
other. There is of course a thermal distribution of such particles above the barrier, and the 
time that it takes the ensemble to lose enough energy to be caught in a well is correspondingly 
distributed. 

The principal ingredients in this calculation are thus 

1. The time that it takes to complete an orbit; 

2. The energy loss per orbit; 

3. The distribution of times at which particles eventually become trapped in one well or 
the other. 

In this section we concentrate on the first two of these ingredients. The third component, 
which then leads directly to the transmission coefficient, is considered in the next section. 
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Because of the slow energy variation it is here convenient to rewrite the dynamical equa- 
tions (|l|) in terms of the displacement and the energy 



(15) 



instead of the momentum O, KM. This simple change of variables leads to 



q = {2[E-Viq)]} 



1/2 



(16) 



E = -27 [E - Viq)] + fit) {2 [E - Viq)]}'/' . (17) 

This set can of course not be solved exactly either, but one can take advantage of the fact 
that for small dissipation 7 the temporal variation of the energy is much slower than that of 
the displacement. 

The first quantity to calculate (approximately) is the time that it takes a particle of 
energy E above the barrier (i? > 1/4) to complete one passage from g = to the edge g+ 
of the potential well and back to g = (cf. Fig. |ip. For this calculation we assume (this is 
the approximation) that the energy E remains fixed during this passage. g+ is given by the 
solution oi E = V{q), i.e., by g+ = (1 + y/AEy^'^. We label t{E) this time of travel between 
q = and q = q+ and back to g = and call this a "half-orbit." The time t{E) is then 
simply obtained by integrating Eq. (p^ : 



m-^r ^ 

Jo {2[E-Vm 



1/2- 



With the potential (0) this can be expressed in terms of a complete elliptic integral of the 
first kind, lC{z): 

^ ' {AEy/^ y aVe j ^ ' 

For low temperatures the particles that determine the behavior of the transmission coefficient 
are primarily those just above the top of the barrier. For these energies an excellent approx- 
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imation to Eq. ( [T9| ) shows that the time to traverse the potential depends logarithmically on 
the excess energy of the particle above the barrier: 

m^\n-^ (20) 
- j) 

It is convenient for this and subsequent results to introduce the notation 

e = E- 1/4, (21) 

that is, the energy above the barrier. In terms of this energy 

16 

te = t{E) ~ In — (22) 
(cf. Fig. 0). The corrections to Eq. ( p2D are of 0{e\ne). 



Equation {2i) is one of the essential results of this section. 

The energy of the particle of course does not in fact remain constant as the particle 
orbits around. Kramers calculated this energy loss in his theory of energy-diffusion limited 
reactions in terms of the action of the particle. We recover his result using the following 
physically transparent argument. To estimate the actual energy loss, say, along a half orbit, 
in the simplest approximation one neglects the fluctuations in Eq. (^) (since particles at 
the high end of the energy distribution primarily lose energy - this overestimates the energy 
loss) and integrates over a time interval t{E). Note that ignoring the fluctuations is not 
tantamount to "ignoring the thermal effects," since the principal thermal effects lie in the 
initial distribution of energies, which we take into account exactly. Subsequently we also take 
account of the fluctuations ignored in this simplest approximation. The q{t) dependence in 
Eq. ([T7|) still poses a problem, though, because it is in general an incomplete elliptic integral 
itself amenable only to numerical integration. However, since q changes rapidly compared to 
E, it would appear reasonable to perform an average of the g-dependent term over a passage 
from g = to This leads to the approximation 

E = -2-fF{E) (23) 
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where 



F{E)^ r dqPE{q)[E-V{q)] (24) 







and Pe{<1) is the probabihty density that q{t) = q. The assumption that this probabihty 
density is in turn proportional to the inverse of the velocity at q(t), i.e., Pe{(i) oc 1/g |T^, 
leads with Eq. (|l^) and appropriate normalization to 

\/2 1 

Pj^(q) = -1— _. (25) 

The integral in Eq. ( P^ can then be expressed in terms of complete elliptic integrals. For the 
energies of interest here an excellent approximation is obtained by retaining the ^'-dependent 
contribution t{E) in Eq. ( pSf ) (which is singular as, E ^ 1/4) but to set = 1/4, the barrier 
height, elsewhere in the integrand. We thus obtain 

Thus, for a first estimate of the energy loss per half orbit one might use Eq. ( pSj) with (pG)): 

47 1 



The next question then is: on the basis of the approximate equation (p7|), how much 
energy does a particle that begins with energy E actually lose per half orbit? The answer 
of course depends on the energy - in particular, according to Eq. (^), it depends on the 
time that it takes that particle to complete its half orbit. Let yu(-E) be the energy loss in a 
half orbit of a particle that has initial energy E. The deterministic equation Eq. ( P?]) can be 
integrated directly and to give in terms of the energy e above the barrier 

— - /iVn — = (e - /i) In + /i. (28) 

3 / £ e 

The solution //(e) of this equation clearly depends on the initial energy e of the particle 
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above the barrier. It is in fact not difficult to see that increases with e. However, this 
dependence is sufficiently mild that the error made in ignoring this variation is no larger than 
the effect of the corrections of 0{e Ine) that have been ignored in Eq. (0). We can therefore 
pick a convenient energy above the barrier to calculate this loss in a half orbit. An upper 
bound to the solution is obtained by setting e >> /i in Eq. (|28|) , which immediately yields 
the simple relation between the energy loss per half orbit and the dissipation parameter 

= ^7- (29) 

This result agrees with that obtained via weak-collision arguments or small dissipation argu- 
ments by a variety of essentially equivalent routes. At the other extreme, we can set e = fi. 
Particles with this energy have just enough energy to complete a half orbit from g = to g+ 
and back to the origin. Solution of Eq. (pSD in this case leads to a /i that is somewhat smaller 
than Eq. (p9D, by approximately 10-15% for dissipation parameters in the energy- diffusion 
limited regime. 

We shall call the value Eq. (^) of fi the "bare" energy loss per half orbit and keep in 
mind that it is an upper bound for fi. 

Equation (WW is the second important result of this section. 



In the absence of fluctuations, the only approximations (to order elne) have been to 
average over a half orbit in calculating the time for the return of a particle to the origin, 
and to assume that in the calculation of this period the energy loss of the particle over the 
half orbit is negligible. Thus, we have approximated the evolution of the energy of a particle 
as determined by Eqs. (p!6|) and (0) with f{t) = by simply assuming an energy loss of 
/i at the end of each half orbit. In Fig. ^ we show the results of the exact integration of 
Eqs. (|16|) and (0) in the absence of fluctuations for a particle starting at g = with an 
energy e = 0.04 above the barrier and moving toward the right. The bottom panel shows 
the trajectory of the particle, q{t). The particle orbits back and forth above the barrier 
several times before getting trapped in the left well. The top panel shows the associated 
phase space portrait. The middle panel shows the decay of the energy (solid curve). The 
decay is not entirely uniform because the rate of energy loss when the particle has high 
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kinetic energy (i.e. as it goes over the deepest portions of the wells) is greater than when 
its kinetic energy is low. However, this non- uniformity is very mild. Our approximation is 
given by the circular symbols: we assume an energy loss of fi per half orbit, and the symbols 
are placed at the times of completion of each half orbit as predicted by our theory. Once 
the particle gets trapped, its orbits become smaller and the energy loss per orbit smaller, 
but these dynamics are in any case not relevant to our problem. We conclude that in the 
absence of fluctuations the assumption of an equal loss of energy fi in the course of each half 
orbit during the recrossing process is appropriate. The other theoretical information used 
in drawing the symbols in the middle panel of the figure are the times it takes to complete 
each of the half orbits. The lower panel in Fig. ^ shows crossings at the times given in the 
second column of Table |I|. Those predicted by Eq. ( P2| ) are presented in the fourth column. 
To obtain these numbers we subtract an energy fi = Aj/S after each half orbit and calculate 
the appropriate for the next half orbit. The agreement is clearly quite good. 

Consider now the evolution of particles in the presence of thermal effects. In Fig. |^ we 
again show the results of the exact integration of Eqs. (|l^) and (p!?!), but now in the presence 
of fluctuations corresponding to a temperature ksT = 0.025. The second and third panels 
represent ensemble averages over 500 particles all starting at g = and moving toward 
the right with an initial energy e = 0.04 above the barrier. The ensemble average is over 
different realizations of the fluctuations. The bottom panel shows the average trajectory 
of the particles, {q{t)). On average, the particles orbit back and forth above the barrier 
several times. About half of the particles eventually get trapped in the left well and the 
other half in the right, whence the average trajectory settles down to values near zero. The 
top panel shows the associated phase space portrait for a typical particle of the ensemble. 
The middle panel shows the decay of the average energy (solid curve). The decay is a bit 
noisier than that of the deterministic case of Fig. ^ because the stochastic ensemble is not 
very large (the fluctuations are of order 0{N^^^'^) where is the number of particles). Our 
approximation (which neglects thermal effects) is again given by the circular symbols: we 
have again assumed an energy loss of fi per half orbit, and the symbols are placed at the 
times of completion of each half orbit as predicted by our theory. Again, the agreement is 
excellent. We conclude that on these time scales even in the presence of fluctuations the 
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Figure 4: Simulation of the evolution of a single particle according to Eqs. (|T^) and (|T^ with 
no fluctuations. The particle starts at g = moving toward the right and with initial energy 
Ei = 0.254 {si = 0.04). The dissipation parameter is 7 = 0.005. Top panel: phase space 
portrait. Middle panel: the decay of the energy (solid curve). Bottom panel: the particle 
trajectory. The times at which the particle crosses g = have been labeled consecutively - 
in this case there are five such times before the particle becomes trapped. Circular symbols: 
the decay of the energy according to our approximation. The symbols are placed at the 
theoretical barrier recrossing times. 

assumption of an equal loss of energy /i in the course of each half orbit during the recrossing 
process is appropriate, and that this value of fi is essentially that obtained disregarding the 
fluctuations. Again, the other theoretical information used in drawing the symbols in the 
figure are the times it takes to complete each of the half orbits. The crossings observed in 
the lower panel in Fig. ^ occur at the times shown in the third column of Table |l|. Note 
that the crossing times are now just a little shorter than in the deterministic simulation, 
an indication that the particles have a bit more energy as they orbit and are thus moving 
a bit faster. This in turn indicates that the energy loss per orbit here is actually a bit 
smaller than in the deterministic case (just barely perceptibly so on the energy scale of the 
middle panels in the figures). We interpret this as the effect of thermal fluctuations that 
counteract the energy loss. On these short time scales the thermal fluctuations barely begin 
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to have an effect on the evolving particles, and hence their behavior is very well captured by 
the dissipative component. We will see subsequently that although the short time behavior 
(such as shown in these figures) is well captured by an energy loss rate that is independent 
of the fluctuations, at longer times the fluctuations play a more important role and the 
corresponding fi should be appropriately renormalized. 




-2.0 I ' ' ' ' ' ' ' ' ' ^ 

0.0 10.0 20.0 30.0 40.0 50.0 

time 



Figure 5: Simulation of the evolution of 500 particles according to Eqs. (0) and ( [LTD with 
fluctuations corresponding to a temperature fcsT = 0.025. The initial energy of each particle 
above the barrier is Ei = 0.04, and the dissipation parameter is 7 = 0.005. Top panel: phase 
space portrait of a typical particle. Middle panel: the decay of the average energy (solid 
curve). Bottom panel: the average particle trajectory. The times at which the particles on 
average cross g = have been labeled consecutively - in this case there are five such times 
before the particle becomes trapped. Circular symbols: the decay of the energy Circular 
symbols: the decay of the energy according to our approximation. The symbols are placed 
at the theoretical barrier recrossing times. 

To estimate the correction to Eq. ( ^91) due to the thermal fluctuations we return to the 
full set Eqs. (0) and (|1^ and consider the Fokker-Planck equation associated with this set 
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for the probability density W{q, E,t\qo = 0, Eq): 



^W = -^^{2[E- Viq)]}'^' W + 27^ [E - Viq)] W 



+ ^ksT^ {2[E - V{q)]Y'" ^ {2[^ - V{q)]Y'" W. (30) 

To find an equation for the evolution of the average energy {E{t)) we multiply Eq. ( pO]) by 
E on the left and integrate over q and E. Upon integration by parts all boundary terms 
vanish and one is left with 

±{E) = -2^{E-V{q))+^kBT (31) 

where 



00 







(/(g, E)) = dE dqW{q, E, t\qo = 0, Eo)f{q, E) (32) 



and the integral over q is over the region where V(q) < E. A well-known generalization of 
the argument surrounding Eq. (|2^) in the energy- diffusion limited regime is to assume that 
the probability density can be separated in the form [IT^ III 



W{q, E, t\qo = 0, Eo) ^ w{E, t\Eo)PE{q) (33) 

(the steps to obtain Eq. ( pH] ) are consistent with this assumption in the special case that w 
evolves deterministically) . The average of Eq. (^) with respect to q is then equivalent to 
the earlier average over a half-orbit at fixed E. With the assumption that we can replace 
{t{E)^^) with t{{E))^^ this finally leaves us with the evolution equation for the mean energy 
that modifies (^) in that it takes thermal fluctuations into account: 

The solution of this equation again clearly depends on the initial energy of the particle, 



20 





Simulation 


Simulation 


Theory 




Deterministic 


ksT = 0.025 




^0 


0.00 


0.00 


0.00 


tl 


6.03 


6.04 


5.99 




12.32 


12.20 


12.16 


h 


18.94 


18.06 


18.56 


u 


26.03 


24.15 


25.24 


h 


34.04 


29.70 


32.33 



Table 1: Zero crossing times of the orbits shown in Figs. § and |. Second column: simulation 
results of Fig. ^ Third column: simulation results of Fig. ^. Fourth column: theoretical 
results. 



and the arguments surrounding the choice (p9D still apply. An exact explicit solution of 
Eq. ( P^D (even to the orders of approximation considered here) is not possible, but a rough 
approximation for low temperatures is 

fi{T) = ^^--fkBT In-. (35) 
3 7 



Of course, numerical integration of Eq. (|3^) is trivial, and in our further analysis we use the 
outcome of such integration. For example, for 7 = 0.005, a value used in a number of our 
simulations, Eq. (|3|) leads to /i(T) = 0.004907 for keT = 0.025 and to fi{T) = 0.004485 for 
ksT = 0.035. These are to be compared with the estimate in Eq. (P^, which is /i = 0.006667. 
Although at short times we have seen that the "bare" value of fi yields excellent agreement 
with simulations, these "renormalized" values do play a more important role at long times, 
as will be seen in the next section. 

The expression or, more precisely, the result of integrating Eq. is the third 

essential result of this section. 
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6 Energy-Diffusion Limited Regime - Transmission Co- 
efficient 



The calculation of the time- dependent transmission coefficient, the goal of this work, is based 
on the principal results of the previous section, namely, the time t{E) = it takes a particle 
of energy e above the barrier to complete a half orbit, and the energy loss /i in the course of 
this time. In order to carry out this calculation, we must keep track of where the particles 
that start above the barrier are at any time before they get trapped in one or the other of 
the wells, and at what time each particle becomes trapped in one well. These quantities of 
course depend on the initial energy of the particle, which is in turn distributed thermally. 

The accounting that leads to the result for the transmission coefficient n{t) turns out 
to be more transparent if we note that due to the symmetry of the problem the following 
interpretation of the terms in Eq. ([I^) leads precisely to the same result as the description 
presented there. We can imagine that all particles are forced to start with a positive velocity 
vq > 0. Then is, as before, the fraction of particles that at time t are over or in the 

right hand well, but now can be interpreted as the fraction of particles that at time t 

are over or in the left hand well. Furthermore, with this interpretation it is clear that 

K_(t) = 1 - (36) 

so it is sufficient to follow one or the other. 

With this interpretation we thus imagine that all particles above the barrier start at 
go = and with a positive velocity Vq distributed according to the Boltzmann distribution. 
Since the energy loss per half orbit is /i essentially independently of the energy of the particles 
(subject to the approximations and conditions discussed in the last section), keeping track of 
particles is conveniently done in "layers" of energy thickness fi. Although the energy loss per 
half orbit is essentially independent of energy, the time for a half orbit is not, since particles 
of higher energy return to g = faster than those of lower energy. This time dispersion must 
be accounted for carefully in constructing K(t). 

Initially the particles are distributed according to a Boltzmann distribution at tempera- 
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ture T. The fraction of particles that are above the barrier is 

P{E > 1/4) = -i- / 6-^/^=^^ = e-^/^'^-^. (37) 

We are interested in keeping track of only this fraction and, in particular, of the fraction of 
this fraction that is over one well or the other as time proceeds. Therefore, we normalize 
this initial fraction to unity, i.e., the distribution of interest is 

p{e)^-^e-'/'-'', 0<£<oo. (38) 
Kb -I 

Initially all particles are over the right hand well, i.e., 

POO 

K+{0) = / dep{e) = 1, ft;_(0) = (39) 

and therefore 



k{0) = 1. (40) 

All those particles that start with an energy between and /i above the barrier never leave 
the right hand well since by the time they return to g = they will have lost energy and 
will thus be below the top of the barrier and unable to move out of the well. The remainder, 
those of energy greater than fi above the barrier, e > /i, do come back and recross the barrier 
to the other side. However, they do so at different times depending on their energy. We can 
write this as follows: 

K4t) = l-^J dse-^/''^^e{t-Q + ---, (41) 

where the Heaviside theta function 6{x) = 1 for a; > and = for a; < 0. The theta function 
insures that recrossing occurs at the right time for that energy and no sooner. The • • • 
indicate that the calculation is not yet complete, since we need to continue to keep track of 
particles that recross back from the left well to the right, and so on. The theta function in 
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Eq. (^) can be implemented by adjusting the limits of integration: 

nJt) = 1 - — - / dee"'^''^^ + ■ ■ ■ , < t < In — , 



1 - -— / dee-'/''^^ + ■■■ , In — < t < ■ ■ ■ . (42) 



Here fi{t) is the inverse of the half orbit time, that is, e = fi{t) is the solution of the relation 

16 

t = t, = In-, (43) 

e 



which immediately leads to 



hit) = 16e-\ (44) 



Note that the upper limit of the range of the time variable in the first line of Eq. ( ^21) is 
precisely the time it takes a particle of energy fi above the barrier to complete a half orbit. 



that is, Eq. (^31) with e = jj,. During this time all the particles with initial energies above 
e = fi will have recrossed the barrier at least once, and all the particles with energies below 
£ = yU are trapped in the right hand well. The upper limits of integration are written as 16 
(rather than oo) because our approximation for (cf. Eqs. ([T9|) and (^OD ) breaks down for 
higher energies, but 16 is, for all practical purposes, infinite when the energy loss per orbit 
is of order 7 and 7 is very small. In our subsequent discussion of the time dependence of the 
transmission coefficient 7 is of order 10~^. On the other hand, the time intervals that define 
the ranges of behavior do involve precisely ln(16//i) because for small fi this is the leading 
contribution to the elliptic integral discussed in Section ^. The integrals in Eq. ( ^2]) can be 
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done explicitly. Treating the upper limit as oo we find 

1 fi 

K+{t) = 1 - e-^^W/'^^^ + ■ ■ ■ , < t < In -, 

^^ 

= 1 - e"^/'^^^ + • • • , In — < t < . . . . (45) 

As before, the ■ ■ ■ indicate that the calculation still continues. 

The next stage of the calculation is to follow the particles that recrossed the barrier once 
and find themselves over the left well. Those that started at t = with energy in the range 
fi < e < 2fi above the barrier and lost energy /i in their first orbit have entered the region 
above the left well with energy between and /i above the barrier, and will therefore never 
make it out of the left well. Only those that started with energy greater than 2/x above the 
barrier and thus recrossed the barrier with energy greater than 1/4 + /i will be able to orbit 
over the left well and return to g = and eventually recross once again to the right hand 
well. The following continuation of Eq. (|4l|) captures this behavior: 

2 poo I roo 

K+{t) = l--— dse-'^'^^^eit + dee-'/'^-^eit - - + ■ ■ ■ . (46) 

Again, the theta functions insure the correct timing of the recrossings and they can be 
implemented by adjusting the limits of integration: 

^^(t) = l--^ / dee-'/''^^ + / &e-^/^^^ + ■ ■ ■ , < t < In — , 



^ dee-'/^-^ + r dee-'/^-^ + ■ • • , In - < t < 2 In ■ 



ksTj^ ksTjj^^,) ' /X (2!)V2;,' 

= 1-^ r dee~''^''^ + ^ r dee-''^^'^ + 21n <t<--- 

ksT J, + ksT A, + , ^ (2!) < * < • 

(47) 

The function f2{t) is the inverse of the sum of the two half orbit times, the first half orbit 
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being that over the right hand well with energy e, and the second that over the left well with 
energy e — fj,. In other words, e = f2{t) is the solution of the relation 

1 16 , 16 

t = t, + t,_^ = In — + In . (48) 

e e — fi 

Exponentiation leads to 

16V* = (49) 
This is a quadratic equation that can easily be solved: 



The upper limit of the second range of the time variable in Eq. (^) is the time it takes a 
particle of initial energy 2^ above the barrier to complete two half orbits, that is, Eq. (^) 
with e = 2^. These particles are trapped in the left well during their second half orbit. The 
integrals can again be carried out explicitly to obtain: 

1 

K+(t) = 1 - e-^^W/'^^^ + e-^^W/'^^^ + • ■ ■ , < t < In — , 



I _ e-f^/ksT ^ ^-Mtyk^T ^ _ _ ln-<t<21n. 



/i (2!)V2^' 



1 

1 _ e-^^/kBT ^ g-2M/A.,T ^ . . . ^ 2 In — ^ < t < • ■ ■ . (51) 



(2!)i/2^ 



The pattern for K+(t) and hence for K,(t) as required in Eq. (|T4D has thus been established. 
We can write the following complete result for the transmission coefficient as a function of 
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time: 



oo 



1 fi 

= 1 + 2 ^(-l)-e-^"W/'=«^, < t < In 

n=l 



oo 



1 _ 2e-^/'=^^ + 2 ^(-l)"e-^"W/^^^, In — < t < 2 In 



n=2 



/i (2!)V2^' 



1 _ 26-'^/'=^^ + 2e-2'^/^^^ + 2 ^(-l)"e-^"W/'=^^, 2 In ^^^^ < t < 3 In 

n=3 



(2!)V2^ (31)1/3;,' 

(52) 

and so on. The function e = fm(t) is the solution of the relation t = + t^-^ + ■ ■ ■ + 
which upon exponentiation turns into the m*'' order polynomial equation (a clear 
generalization of the results Eq. (H3) for m = 1 and Eq. (HO) for m = 2): 



(16)™e-* = e[e- ij]---[e-{m- l)/x]. (53) 

The solution of Eq. (^3|) can in general not be found in closed form for m > 3. However, an 
excellent approximation is 

fm{t) ^[m- (m!)^/'"] fi + lee-*/'". (54) 

This form is exact for m = 1, and it is exactly correct for all m at the upper limit of the 
time range that defines the trapping of the particles whose initial energy is e = m/x. In other 
words, fm(t) as given in Eq. ( ^4]) is exactly correct at the particular time t = t^+t2fj,+- ■ ■+tmfi- 
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Equation ( |5^ and its asymptotic limit 



oo 



(t ^ oo) = l + 2^(-l)"e 



nii/ksT 



n=l 



tanh 



( 



2kBT 



) 



(55) 



are the principal results of this paper. These results generalize those of KT in the Markovian 
regime to the energy-diffusion limited case. 

It is useful to examine the information contained in the various contributions to Eq. (^) 
and the preceding pieces that were used to construct it. The first line contains all first 
recrossings of the barrier by all particles that start with sufficient energy above the barrier 
to recross is at least once {e > fi). The function /i(t) accounts for the fact that this first 
recrossing occurs at different times for particles of different initial energy, the last ones (those 
closest to the barrier that do make it around) recrossing at time t = ln(16//i). The particles 
of higher energy recross first because they are orbiting at higher velocity, and may be back 
for their second recrossing while those of lower energy are still awaiting their first crossing. 
The time distribution of the second recrossing is contained in f2{t)- Even while this is 
occurring, some particles might already be undergoing their third recrossing, as contained 
in fslt). There are fewer and fewer of these faster particles because of their initial thermal 
distribution - this information, too, is contained in the exponential factors. Of course while 
all these recrossings are going on, the particles are losing energy and, depending on their 
initial energy and how often they have gone around, they become trapped at sequential times 
as expressed in the subsequent hues of Eq. ( p^ ) . 

We note the nonlinear dependence of on the energy loss parameter fi and hence on 
the dissipation parameter 7. The Kramers and other weak collision results in the energy- 
diffusion limited regime correspond to the retention of only the leading term in a small-7 
expansion of our result, which then yields ~ •y/kBT. 

The detailed behavior of the time-dependent transmission coefficient as given by Eq. (|5^ , 
and its asymptotic value Eq. (|5BD, are compared with simulations in the next section. 
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7 Discussion of Results and Comparisons with Numer- 
ical Simulations 



Figures § and show the transmission coefficient K,{t) as a function of time for two sets 
of parameters. In Fig. ^ the dissipation parameter is 7 = 0.005 and the temperature is 
ksT = 0.025; in Fig. |^ they are 7 = 0.01 and ksT = 0.05. Note that these are the only 
input parameters. The sohd curves are simulations of an ensemble of 4000 particles. Each 
particle follows the equations of motion (|l]), and the transmission coefficient is calculated 
from the resulting trajectories according to the reactive flux formalism described earlier. 
The distinctive features of the time dependence are 1) the time at which the transmission 
coefficient drops rather abruptly from its initial value of unity, and the slope of this drop; 2) 
the frequencies and amplitudes of the oscillations; and 3) the asymptotic value, identified as 
the equilibrium value Hgt of the transmission coefficient. The dashed curves are the result 
of our theory, Eq. (|52D, with the bare value /i = 47/8 for the energy loss per half orbit. We 
stress that there are no adjustable parameters in these curves. The agreement between the 
theory and simulations is clearly very good in both cases. Without adjustable parameters, 
the theory captures each of the three distinctive features listed above. Indeed, we stress that 
with a single expression we are able to reproduce the temporal behavior of the system over 
essentially all time scales. 

First we discuss those aspects of the theory that work extremely well, and then we suggest 
ways in which improvements could be made if even closer agreement is desired. The theory 
essentially exactly captures the ffist striking feature in the transmission coefficient, namely, 
the time and slope of the ffist abrupt drop. The drop corresponds precisely to the point at 
which the KT theory and therefore also the Grote-Hynes theory deviate from the simulation 
results in Fig. (the parameter values there and here are different). The time of the drop 
and the slope of the ffist drop are determined by the particles that ffist recross from one 
well to the other. Although this time behavior is contained in the full sum of the ffist line of 
Eq. (I52D , the ffist portion of the sum, that is, the terms shown explicitly in Eq. (|45|) , already 
capture this behavior very well. These terms account for the first crossing over the barrier 
and do not include the effects of recrossings (which are contained in the higher fm{t), as 
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Figure 6: Transmission coefficient of an ensemble of particles at temperature ksT = 0.025 
and dissipation parameter 7 = 0.005. The particles are initially distributed above the barrier 



according to Eq. (|37|). Solid curve: simulation of 4000 particles. Dashed curve: our theory, 
Eq. (p2D, with the bare value /i = 47/3 = 0.006667. Thick line that intersects the right 
vertical axis: value of the equilibrium transmission coefficient Kgt obtained from Eq. (^) 
using the value of 11 renormalized by thermal fluctuations, which for this temperature is 
fi{T) = 0.004907. 



described earlier). In other words, the drop is principally due to the particles near the top 
of the barrier that manage to recross. The determining features in the early time behavior 
of K,{t) are thus the time range < t < ln(16//i), the value of fi, and the decay function 
fi{t). Once the least energetic particles are trapped, however, the particles of initially higher 
energies (although fewer in number) sequentially dominate the oscillatory behavior of /t(t), 
that is, the higher terms in the sum become sequentially important in the time dependence 
of K,{t) as time increases. The subsequent oscillations thus reflect recrossings of the barrier. 
These are also captured quite accurately. The simulations and the theory settle to asymptotic 
values that are extremely close to one another. 

The theory could of course be improved in various ways. First, the theoretical oscillations 
are sharper and continue for a longer time than those of the simulations because the theory 
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Figure 7: Transmission coefficient of an ensemble of particles at temperature = 0.05 
and dissipation parameter 7 = 0.01. The particles are initially distributed above the barrier 



according to Eq. (|37|). Solid curve: simulation of 4000 particles. Dashed curve: our theory, 
Eq. (|52D, with the bare value fi = 47/3 = 0.01333. Thick line that intersects the right 
vertical axis: value of the equilibrium transmission coefficient Kst obtained from Eq. (^) 
using the value of /j, renormalized by thermal fluctuations, which for this temperature is 
/i(T) = 0.007892. 



averages or altogether ignores additional effects that lead to particle dephasing. For example, 
the sharp (deterministic) time cuts (which cause the sharp minima and maxima in the theory) 
are in reality somewhat distributed due to the fact that the particles lose energy continuously 
and not just at the ends of half orbits as assumed in the theory, and due to fluctuations that 
allow energy gains throughout the dynamics. Indeed, particles that we take to be trapped 
forever can in reality recross the barrier due to thermal fluctuations (but in very small 
numbers at low temperatures, as reflected in the small value of the equilibrium transmission 
coefficient). Furthermore, the average energy loss per half orbit, fi{e), has been assumed 
independent of e, but in reality fi increases with e, even in the "bare" approximation that 
ignores thermal fluctuations in the calculation of fi. This dependence is contained in Eq. (pSf), 
but we have approximated the solution of Eq. (pSf) by the e- independent Eq. (^9]), which is 
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an upper bound on /i. On the other hand, thermal fluctuations contribute to a decrease in 
/i as a function of temperature. Indeed, if one analyzes average energy loss trajectories such 
as the one shown in the middle panel of Fig. |^ but for slightly higher initial energies, one 
finds that there is a small upward curvature that indicates that fi starts with the bare value 
but that it becomes renormalized by thermal fluctuations as time proceeds. 

The theoretical equilibrium values of the transmission coefficient achieved by the dashed 
curves in Figs. |^ and |^ are a little higher than the simulation result. Note that the entire 
dashed curve in both cases has been calculated using the upper-bound bare energy loss 
parameter /i = 47/3. If we rely on our result Eq. (|55D , the small difference between simulated 
and calculated results would indicate that the actual value of fi at long times is a bit lower 
than the value assumed in the dashed curve. We also show in the figures the equilibrium 
transmission coefficient calculated in each case according to Eq. (|55D , but now with the 
renormalized energy loss parameter that follows from solving Eq. (|34D. For the parameters 
in Fig. I this value is i2{T) = 0.004907, and in Fig. it is /i(r) - 0.007892 (thick short lines 
that intersect the right vertical axis). These values of the stationary transmission coefficient 
agree extremely with the simulation results in both cases. 

Finally, it is interesting to pursue in more detail the dependences of the equilibrium 
transmission coefficient Hgt on the input parameters (7 and ksT) and, in particular, to 
evaluate the prediction Eq. (|55| ) for these dependences. As noted earlier, the existing results 
correspond to a small-7 expansion of our result. In practically every instance that we have 
found in the literature, the equilibrium transmission coefficient vs dissipation parameter is 
drawn on a log-log plot that would mask the dependence that we are focusing on. For this 
analysis we return to Fig. |^ and, in particular, the low dissipation portion of the figure. We 
present two sets of results for vs 7 on a semi-log plot, one for temperature ksT = 0.025 
(the temperature used throughout most of this paper) and the other for ksT = 0.05. The 
solid circles are the simulation results for the lower temperature, and the triangles are for 
the higher temperature. The first and very clear point to note is that the 7 dependence 
is indeed not linear beyond the very lowest dissipation parameters. The upward curvature 
with increasing 7 is clear. We then ask whether the tanh dependence in Eq. (E^) captures 
this curvature. The only small uncertainty lies in the value of the parameter /i to use 
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in answering this question. For the particular dissipation parameter value 7 = 0.005 and 
temperature ksT = 0.025, and also for 7 = 0.01 and fc^T = 0.05, we have established that 
the equilibrium value we predict is quite accurate if we use the bare fi and even more accurate 
if we use the renormalized value fi{T). In Fig. |^ we show the theoretical curve Eq. (55) for 
two temperatures using the renormalized value fi{T). The solid curve is for temperature 
ksT = 0.025, and the dashed curve for fc^T = 0.05. Clearly the theoretical curves capture 
the simulation results very well over a substantial range of dissipation parameters values. 



8 Conclusions 

We have developed a theory for the time dependent transmission coefficient in the energy- 
diffusion limited regime of the Kramers problem. Our work complements that of Kohen and 
Tanner whose theory covers only the diffusion limited regime. We arrive at an explicit 
analytic prediction for the transmission coefficient that describes the rate at which particles 
in the potential of Fig. |l] move over the barrier from the region associated with one well to 
that associated with the other. Our result depends only on the two input parameters of the 
problem, the dissipation and the temperature, and involves no adjustable parameters. For 
low dissipation and low temperature, our result captures the correct behavior on all time 
scales, leading to excellent agreement with simulations in the short time regime where the 
transmission coefficient oscillates, and also in the long time regime where it settles to an 
asymptotic equilibrium value. Our prediction for the equilibrium value as a function of the 
dissipation parameter and temperature extends beyond the regime of existing theories in 
this regime. 

The theory starts with the Langevin equation for the system and proceeds in two stages. 
First, we calculate trajectory features such as the time it takes a particle to complete a half 
orbit and the energy loss in the course of a half orbit. This provides the ingredients that 
allow us to determine when particles cross and recross the barrier between one well and the 
other as a function of their initial energy, and how long it takes for a particle to become 
trapped in one well or the other. The second stage is the combination of these ingredients for 
an ensemble of particles that are initially distributed canonically. The resulting oscillatory 
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early time behavior of the transmission coefficient reflects barrier recrossing events of initially 
energetic particles. The equilibrium value of the transmission coefficient reflects the residual 
steady state flux of particles essentially trapped in one well or the other. 

We also note that we have carried out new simulations in regimes where theoretical pre- 
dictions were available but had not been checked, and in the regimes of interest in this paper 
where neither theory nor simulations were previously available. Through these simulations 
it has been possible for us to check our own work as well as that of others (speciflcally, the 
KT results for the diffusion-limited regime), and we now have available simulation results 
that can serve as a backdrop for further theoretical developments. 

That a number of chemical reactions occur in the low dissipation regime has, of course 
long been clear, not only on the basis of theoretical considerations but on the basis of a 



variety of experiments (see e.g. [|l^, |T6[). The ability to probe reactions on a very short time 
scale is much more recent [^, ^. The theory we have presented provides a more complete 
understanding of energy- diffusion limited reactions than has heretofore been available. To 
further complete the picture, our work continues in order to generalize our model, most 



immediately to non-Markovian processes [|I7 
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